Direct measurement of the magnetic penetration depth by magnetic force microscopy 
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We present an experimental approach using magnetic force microscopy for measurements of the 
absolute value of the magnetic penetration depth A in superconductors. A is obtained in a simple and 
robust way without introducing any tip modeling procedure via direct comparison of the Meissner 
response curves for a material of interest to those measured on a reference sample. Using a well- 
characterized Nb film as a reference, we determine the absolute value of A in a Ba(Feo.92Coo.os)2As2 
single crystal and a MgB2 thin film through a comparative experiment. Our apparatus features 
simultaneous loading of multiple samples, and allows straightforward measurement of the absolute 
value of A in superconducting thin film or single crystal samples. 
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I. INTRODUCTION 



The superconducting coherence length (£), magnetic 
penetration depth (A), and their anisotropy are fun- 
damental parameters that characterize superconducting 
materials^— A number of important properties, such as 
superconducting critical fields and superconducting fluc- 
tuations that affect vortex dynamics, can be obtained if 
the parameters above are known* 1 ^ The value of £, which 
depends on the Fermi velocity and the condensation en- 
ergy of the superconducting state, can be estimated from 
the upper critical field (i? C 2) using the Ginzburg-Landau 
theory. A is related to the density of superconducting 
electrons,— and, in contrast to £, precise determination 
of its absolute value is notoriously difficult 1 -^ due to de- 
magnetizing effects, topography-related surface barrier, 
and inhomogencity of the sample. Typically A is calcu- 
lated by indirect methods. Several experimental tech- 
niques such as tunnel diode oscillator (TDO)j£ temper- 
ature dependence of the flux explusion in the Meissner 
stated or of the reversible magnetization in the mixed 
state using superconducting quantum interference device 
(SQUID) magnetometers j 9 - mutual inductance^ surface 
impedance fii infrared reflectivity^ muon spin resonance 
(/iSR) ) 13 ' 14 nuclear magnetic resonance (NMR) and 
magnetic force microscopy (MFM) 1 ^— have been em- 
ployed for measurements of A in thin films and bulk sam- 
ples. Each of these methods has its own advantages as 
well as limitations, and some of them require simulations 
with multiple fit parameters. For example, TDO mea- 
surements depend on the quality and thickness of an Al 
film deposited on top of a superconductor, which may 
not be fully reproducible yielding errors in the obtained 
values of A. The temperature dependent Meissner re- 
sponse method is only sensitive for thin films with mag- 
netic filed parallel to the surface. The reversible mag- 
netization method works only in clean samples or ma- 
terials having an extended vortex liquid phase. Mutual 
inductance and /iSR techniques are limited to thin films 



and bulk samples, respectively. /J.SR measures the sec- 
ond moment of the magnetic field distribution around 
vortex, consequently details of the vortex structure and 
the muon's location affect experimental accuracy. In- 
frared reflectivity allows measurement of the anisotropy 
A by polarization of an incident light. 

MFM has been widely used for studies of super- 
conductors, particularly for imaging and manipulation 
of vortices in superconducting thin films and single 
crystals^ - — Recently, MFM was also used as a local 
probe of the magnetic penetration depthj 1 ^— where the 
values of A were extracted by fitting either an MFM signal 
from a single vorte:*i£ (note, this method doesn't work 
for systems with large values of the magneic penetra- 
tion depth 2 ^) or a Meissner response . 18 ' 19 These meth- 
ods require a thorough characterization of the probe tip 
and a well-defined simulation model (dipole-monopolc, 
monopolc-monopole, etc.) describing interactions be- 
tween the MFM tip and a superconducting sample. How- 
ever, modeling procedures with multiple fitting parame- 
ters introduce uncertainties in the resulting A values. In 
spite of technical difficulties MFM has certain advantages 
e.g. localization of measured A values providing a route 
to explore anisotropy of A throughout the sample. In our 
approach, we obtain local values of A by directly com- 
paring the Meissner curves for the sample under inves- 
tigation with those obtained for a reference sample with 
a well-known A; we emphasize that the measurements 
are done using the same MFM tip during the same cool- 
down. We observed strong dependence of the Meissner 
response curves on the shape of the MFM tip. When 
the tip crashes or even slightly touches the sample sur- 
face, Meissner response curves can change significantly. 
During our experimental procedure we verify that the 
MFM tip does not change its properties by comparing 
the Meissner curve from a Nb reference sample before 
and after the measurement. In this paper we demon- 
strate the validity of our method by determining A in 
a Ba(Feo.92Coo.o8)2As2 single crystal and a MgB 2 film. 
Our results demonstrate that the same procedure can be 
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used in any single crystal or thin film superconducting 
samples with a thickness greater than A (for thinner sam- 
ples A can be corrected in a straightforward manner . 26 ' 27 ) 

II. EXPERIMENT 

The measurements described in this paper were 
performed in a home-built low-temperature MFM 
apparatus^ We have developed an additional capabil- 
ity of mounting multiple samples (including a reference 
sample), as shown in Fig. [2(a), for acquiring a complete 
set of MFM data for each of the samples within a sin- 
gle cool-down. The absolute value of A is obtained by a 
simple comparison of the Mcissncr curves for the Nb ref- 
erence to those obtained in the sample of interest. The 
Meissner response curve is first measured in a homoge- 
neous reference sample (Nb film) as a function of the 
tip-sample separation. Then the cantilever is moved over 
the sample of interest, and its Meissner response curve 
is obtained. Direct comparison of these curves (compar- 
ative experiments) yields the absolute value of A in a 
sample under investigation. The value of A in the refer- 
ence sample (Nb film) was verified by a different MFM 
technique and a SQUID magnetometry measurement f^ 7 - 
The reference Nb thin film (T c w 8.8 K, where resistance 
drops to zero) has a thickness of 300 nm and was grown 
by an electron beam deposition. A Ba(Feo.g2Coo.os)2As2 
single crystal (T c w 22 K from specific heat capacity mea- 
surements) was grown out of FeAs flux . 29 ' 30 The 500-nm 
thick MgB2 film (T c « 38.3 K, zero resistance tempera- 
ture) was grown by reactive evaporation* 3 -! The MFM 
measurements were performed using a high resolution 
Nanosensors cantilever— that was polarized along the tip 
axis in a 3 T field of a superconducting magnet. The su- 
perconducting samples are zero field-cooled for the Meiss- 
ner experiment and are field-cooled in a field of few Oer- 
sted, applied perpendicular to the film surface and paral- 
lel to the probe tip, for imaging of vortices. All samples 
are electrically grounded to eliminate a possible electric 
force contribution from a stray charge on the sample to 
the magnetic Meissner force. 

III. RESULTS AND DISCUSSION 

A. MFM measurement in the Nb film 

As the magnetic tip approaches a superconducting 
sample, it experiences a Mcissncr force (Meissner re- 
sponse) induced by the shielding currents in the sample 
that screen the magnetic field of the tip. The experimen- 
tal procedure is as follows: First the probe tip is brought 
close to the reference sample [Nb film, position 1 in panel 
(a) of Fig. [T] and the Meissner response is recorded as a 
function of the tip-sample separation at a certain temper- 
ature T . There should be no vortices present in a large 
field of view of the sample, i.e., the sample should remain 




FIG. 1: (Color online) (a) Sample holder with multiple sam- 
ples. The Nb thin film (300 nm) and the MgB 2 thin film 
(500 nm) samples are labeled as 1 and 2, respectively, (b) 
Schematic illustration of the Meissner experiment. The Meiss- 
ner response force between the probe tip and the sample can 
be regarded as the force between the real tip and image tip 
at 2(z + A), when z > A.— 

in a pure Meissner state. Figure GJa) shows the Meiss- 
ner state with no vortices present in the 25 fim x 25 [im 
field of view after a small compensation field was applied 
above T c of Nb to compensate the remnant stray field of 
the superconducting magnet. 

Figure [2Kb) shows several Meissner response curves for 
the Nb film reference sample between 5.5 K and 8 K. 
The Mcissncr response force is a function of A and the 
tip-sample separation z, Fuj = F{z + A) i 18 ' 33 For a given 
temperature T, F M [z + A(T)] = F M [z + A(4 K)+SX], 
where SX = A(T) — A(4 K). Therefore, to determine the 
A value at a particular temperature T, Fm[z + A(T)] is 
shifted along the z axis to coincide with Fm[z + A(4 K)], 
a reference curve measured at 4 K, and the value of the 
shift yields the value of S\(T). 

The Meissner response curve taken at 8 K shows a be- 
havior very different from the data taken at lower tem- 
peratures. Figure[Uc) shows the MFM image acquired at 
4 K after the Meissner measurement was performed at 8 
K. The scan areas in Fig. [21(a) and Fig.[2](c) are the same. 
The vortices seen in Fig. [2tc) were generated by the field 
of the probe tip. They cause magnetic field leakage and 
weaken the Meissner response (see Fig.[2]Jb), 8 K curve). 
Such a behavior is observed when A is comparable to the 
film thickness, and the magnetic field from the tip can 
not be fully screened] 26 ' 27 Special care must therefore be 
taken with the Meissner technique in the vicinity of T c 
in thin films since increasing A may take the sample out 
of the pure Meissner state. 

To verify the 4 K value of A(4 K)= 110±10 nm for 
the Nb reference, measured previously using a differ- 
ent MFM technique,— we performed measurements of 
A using the SQUID magnetometry. The Nb reference 
sample used in MFM measurements (film L=3.2 mm 
x W=4.2 mm) was oriented carefully with H paral- 
lel to the surface along the side L, using a home-built 
sample holder. By measuring the transverse component 
m± of the magnetic moment m in the Meissner state 
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FIG. 2: (Color online) (a) MFM image of the Nb sample at 
T = 4 K; no vortices are present in a 25 /im x 25 ^im field 
of view, (b) Meissner response curves taken over the scan 
area shown in (a), (c) The MFM image taken with the same 
field of view as in (a) after the Meissner response experiment 
at 8 K; vortices are nucleated by the magnetic field of the 
MFM tip. (d) Temperature dependence of the penetration 
depth A(T) in Nb obtained by overlaying Meissner curves on 
top of the 4 K reference curve (black circles). Green squares 
represent SQUID data; the red-dashed curve is a fit to the 
BCS model. 



(which should be zero in the case of perfect alignment) 
we confirmed that the miss-orientation between H and 
the film surface was <fi m is ~ 0.2°. In this configura- 
tion, in the Meissner state the component of m paral- 
lel to H is mil = (H/A-k) x LWd e ff, where the effective 
thickness d eff (T) = d - 2A(T) tanh(d/2A(T)) is smaller 
than the geometrical thickness due to the field penetra- 
tion from both surfaces i£ We measured my versus H 
at several T, and from the slopes dm^/dH we extracted 
A(T). The main source of error in this method is the 
spurious contribution to mil due to the projection of the 
Meissner signal arising from the transverse field compo- 
nent, ~ m±(f> m i Sl which is ~10 % of my, thus introduc- 
ing an error of ~ 5%. However, it does not distort the 
functional dependence of A(T). The absolute values of 
A(T) from these SQUID measurements are marked with 
green squares in Fig. [SJd) and agree well with both the 
MFM data (black circles) and the isotropic single-gap 
BCS model (red-dashed curve). 
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FIG. 3: (Color online) Meissner response curves obtained 
from (a) the Nb reference, (b) a Ba(Feo.92Coo.os)2As2 sin- 
gle crystal at 4 K. Different slopes of the Meissner curves 
obtained from each sample indicate a systematic change of A. 
Inset: The Meissner curve of Ba(Feo.92Coo.os)2As2 is shifted 
by 2=200 nm along the z axis to overlay the Meissner curve 
of the Nb reference sample. The difference of the penetration 
depths AA= 200 nm can be obtained from the value of the 
shift along the z axis. 



B. Measurements of the absolute values of A in a 
Ba(Feo.92Coo.o8)2 As2 single crystal 

In recent years iron-based pnictide superconductors 
have drawn a great deal of attention since these sys- 
tems exhibit superconducting properties which are in- 
termediate to conventional BCS superconductors and 
high-T c cuprates. Magnetic penetration depth of a 
Ba(Fei_ 2 ,Co 2 ,)2As 2 system (BFCA), the so-called 122 
family, was investigated using variety of techniques in- 
cluding TDOri M SR^ and MFMJ£ The Stanford group 
utilized an MFM technique and reported absolute val- 
ues of A in this system as a function of doping level. 
The authors used fitting algorithms to approximate the 
tip magnetization and calculate the values of \J£ We 
applied our direct (comparative) technique to measure- 
ments of A in the very same system [ Ba(Fei_ a; Co a; )2As2 
(x=0.08)] to demonstrate the validity of our approach. 
Figure [3] shows Meissner curves as a function of the tip- 
sample separation z obtained from BFCA (blue squares) 
and a Nb reference (red circles), respectively. The slow 
decay of a frequency shift in BFCA sample compared 
to that in the Nb reference indicates greater values of 
A in the BFCA sample. The expression for the Meiss- 
ner response force in Nb (assuming monopole-monopolc 
interaction between the tip and the sample) can be writ- 
ten as follows: F^ issner = r^^ss , where A is a pref- 
actor that reflects the sensor's geometry and the mag- 
netic moment, $ is a single magnetic flux quantum, 
z is the tip-sample separation, and A^ b is the in-plane 
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FIG. 4: (Color online) Single vortex images in (a) the Nb thin 
film and (b) the MgB2 thin film samples at T = 4 K. Both 
images were taken with a tip-lift height of 400 nm. (c) The 
single vortex profiles taken along the dotted lines in (a) and 
(b). (d) Meissner response curves taken at 4 K in Nb (black 
square) and MgB2 (red circle) with the same experimental 
condition. The green-crossed marks represent that the red 
circle (MgB2), after shifted by 90 nm along the z axis, is 
overlaid over the black square (Nb) to show the validity of 
our approach. 



magnetic penetration depth since the shielding current 
runs within the basal plane. The Meissner force in 
BFCA has the same functional form but different X BFCA : 
F MeSaner = (z+ xB%A y , - Two Meissner curves become 
identical (F*iP s A ner = F™>. ) when the tip lift z com- 
pensates for differences \ BFCA — \ Nb = SX. In other 
words, by shifting the Meissner curve for BFCA along 
the z axis to overlay the Meissner curve measured for 
Nb, the SX value can be extracted. By adding the shifted 
value Sz = 5X = 200 nm to the X Nb (110 nm) one ob- 
tains A SFCA =310±30 nm. This value is close to the one 
reported previously— The inset in Fig. [3] demonstrates 
that two Meissner curves overlay each other very well 
after shifting the BFCA curve along the z axis by 200 
nm. 



C. Measurement of the absolute value of A in a 
MgB 2 film 

We used the same approach to measure the absolute 
value of A in a MgB2 thin film sample. The Meissner 
response curve taken at 4 K in the Nb was used as a 
reference Meissner curve and compared to the Meissner 
response curve measured in MgB2. The offset between 
these two curves, shown in Fig. 2Jd), yields the absolute 
value of A at 4 K. The SX between the Nb and MgB 2 
curves at 4 K equals 90 nm yielding A(4 K)= 200±30 
nm in MgB2. Our experimental error is 10% — 15%, and 
depends on the magnitude of A and the system noise 
level. Absolute values of A in MgB 2 measured with var- 
ious techniques range from 40 nm to 210 nrni— The 
large variation of A in MgB 2 may be due to inclusion of 
impurities, such as C, N, and Al, which replace either 
Mg or B and significantly affect the electronic structure 
of the system due to the two-band nature of MgB 2 . The 
origin of the impurities-induced large value of A and its 
temperature dependence obtained from the Meissner re- 
sponse curves in the MgB 2 film taken at different tem- 
peratures (not shown) will be described elsewhere^ 3 - 

We also imaged individual vortices in both Nb and 
MgB 2 . Direct comparison of the vortex profiles provides 
additional information on the magnitude of A. The maxi- 
mum force gradient at the center of a vortex, max(<9//<9z) 
(MFM is sensitive to a force gradient), is proportional to 
(z + Xab)^ 1 ^ 3, for a monopole-monopole model of the tip- 
vortex interaction ! 23 ' 24 ' 44 In this model, the larger value 
of A results in a smaller force gradient at the center of 
a vortex. Figures|4|a)-(b) show well-isolated vortices ac- 
quired in the Nb reference and the MgB 2 sample at 4 
K. We estimate the magnetic field to be no more than 
0.2 Oe based on the field calibration. Direct comparison 
of vortex profiles [see Fig. 0Jc)], taken along the dotted 
lines in Figs. Ufa) and (b), shows that the force gradi- 
ent in the MgB 2 sample is smaller than that in the Nb 
sample. This indicates a larger value of A in our MgB 2 
film. It is worth noting that the relatively high tip-sample 
separation (400-nm tip-lift) and the tip geometry are re- 
sponsible for the broadening of the vortex force profiles 
in both Nb and MgB 2 samples in Fig. |U[c). 



IV. CONCLUSION 

In conclusion, we have developed an experimental 
method and apparatus to determine the absolute value 
of the magnetic penetration depth A in superconduct- 
ing samples by comparing their Meissner response curves 
to those acquired for a homogeneous Nb reference film. 
We used this method to obtain the absolute value of A(4 
K)=310 ± 30 nm in a Ba(Feo.g 2 Coo.os)2As 2 single crys- 
tal, consistent with the value reported by the Stanford 
group. We also measured A(4 K)= 200 ± 30 nm in a 
MgB 2 film. The large A comes from the nature of the 
two band superconductivity, and from inclusion of impu- 
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rities such as C and N. Our MFM apparatus allows us 
to simultaneously load and investigate multiple samples 
(over ten samples can be studied at once, providing an 
opportunity to explore the complete phase diagram of a 
superconducting system), and most importantly to use 
the same cantilever tip for both Nb reference and the 
samples under investigation in a single cool-down. This 
capability enables in- situ calibration of the MFM tip on 
a known homogeneous Nb sample and does not intro- 
duce any additional uncertainties due to modeling of the 
tip geometry and the resulting tip field. The validity of 
our approach is established by comparing the MFM and 
SQUID magnetometer measurements of the temperature 
dependence of A in the Nb reference film. Our experi- 
mental approach opens the possibility of measuring the 
absolute value of A(T) in film and bulk superconducting 



samples. 
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